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Abstract 



The mass of the bottom quark and the strong couphng constant a^ are 
determined from QCD moment sum rules for the T system. Two analyses 
are performed using both the pole mass Mb as well as the mass rrn, in 



the MS scheme. In the pole-mass scheme large perturbative corrections 



resulting from coulombic contributions have to be resummed. In the MS 
scheme this can be avoided by an appropriate choice for the renormalization 
scale. For the bottom quark mass we obtain M^ = 4.60 it 0.02 GeV and 
iTibiiT^b) = 4. 13 lb 0.06 GeV. Our combined result from both determinations 
for the strong coupling is as{Mz) = 0.119 it 0.008. 



1 Introduction 

The Upsilon system constitutes a rich source of information about the strong 
interaction dynamics. The bottom quark mass is sufficiently heavy for a non- 
relativistic description to be a good starting point to analyze the quark-antiquark 
forces. Thus, potential models — including relativistic corrections — have been 
successfully used to understand the spectroscopy of the corresponding mesonic 
bound states |l]-§[. At the same time, the small size of the hadronic system makes 
possible to attempt a short-distance approach. At the relevant energy scale, a^ is 
small enough to allow (at least for the lowest levels) a well-grounded quantum field 
theory analysis with perturbative QCD tools, including non-perturbative correc- 
tions through the Operator Product Expansion (OPE) [Q. While the coulombic 
part of the bb potential is obviously related to the static piece of the gluon- 
exchange interaction, a systematic short-distance investigation provides a better 
understanding of the remaining terms in the heavy quark potential, in terms of 
fundamental QCD parameters jBHTS . 



The short-distance description in terms of quarks and gluons is specially well 
suited for inclusive quantities, where no reference to a particular hadronic bound 
state is needed. The vacuum polarization n(g^) induced by the heavy-quark 
vector current 67^6 is then a key ingredient to investigate the J^ = 1^ bb states. 
Its imaginary part can be experimentally determined from the e^e^ -^ bb cross- 
section: 

Rb{s) = QlR{s) = ^^"^ ": ^^^\, = 127rQ^inin(s + ze). (1.1) 

cr(e+e —>■ jj,^ jj, ) 

On the other side, n(g^) can be calculated theoretically within the OPE. To keep 
our equations more general, let us consider the vector current j^(a;) = {Q'~iij.Q){x), 
Q{x) being a heavy quark field of mass M, specifically the bottom quark in our 
case: 

(q,qu-9,.q^)Ti{q') = t Jdxe''^' {T{j,{x) j^}) ■ (1-2) 

Throughout this work, M corresponds to the pole of the perturbatively renor- 



malized propagator, whereas the running quark mass in the MS scheme [[16 
renormalized at a scale /i will be denoted by m{^). 

Using a dispersion relation the nth derivative of n(s) at s = can be expressed 



in terms of the nth integral moment of R{s): 



2 / , \ n oo 







mrids^. (1.3) 



.=0 ^ 



For later convenience, the moments A^„ are defined to be dimensionless quan- 
tities. In addition, it will prove useful to express the moments A^„ in terms of 
integrals over the variable v = Jl — AAf^/s, 

1 

Mn = 2 jdvvil-v^y-^R^v). (1.4) 



Under the assumption of quark-hadron duality, the moments M.n can either be 
calculated theoretically in renormalization group improved perturbation theory, 
including non-perturbative condensate contributions, or can be obtained from 
experiment. In this way, hadronic quantities like masses and decay widths get 
related to the QCD parameters a^, m^ and condensates. 

For large values of n, the moments become dominated by the threshold region. 
Therefore, they are very sensitive to the heavy quark mass. This fact has been 
exploited since the very first QCD analyses of charmonium and bottomium [^ 
p!7|-p3| to extract rather accurate values of iric and mi,. More recently, it has been 



suggested by Voloshin ||2J] that the large-n moments can also be used to get a 
precise determination of a^ from the existing data on T resonances. 

The perturbative calculation of the moments contains powers of as\/n [|17|,|18 



25|] , which correspond to the coulombic contributions; they are associated with 
the near-threshold quark-antiquark configurations at typical velocity v ~ l/v^) 
so that as\fn ~ asjv is the familiar Coulomb parameter. At large n these 
coulombic {as\fnf- terms should be explicitly summed up to assure a reasonable 
convergence of the perturbative series. By the same token, this large-n behaviour 



implies a big sensitivity to the value of a^ |]2^ 



In ref. |^J] the large-n moments A^„ have been studied with a non-relativistic 
expansion in powers of \jn. Fitting the OiXjn) contribution from the sum rules, 
the analysis of the moments n = 8, 12, 16 and 20 gave the result: M^ = 4827 ± 
7MeV and o^'^{Mz) = 0.109 ± 0.001. The quoted errors are claimed to include 
the experimental uncertainties and the theoretical uncertainty due to subleading 



1/n terms p4| 



The reasoning of ref. [Q looks indeed very suggestive. The large-n moments 



are dominated by the first T resonance. Thus, one is actually starting with 



a confined bound state. In spite of tliat, our ability to make an explicit sum 
of the coulombic contributions allows to make an impressive determination of 
the perturbative coupling. Obviously, an accurate analysis of the theoretical 
uncertainties is called for. 



The analysis of ref. [|J] was performed at 0{as), i.e., only the 0{1) and 0{as) 
perturbative contributions to the correlator n(g^) were included. Therefore, the 
scale and scheme dependence of as was not under control. Adopting the MS 
scheme, the running of «« was included in the Coulomb potential, and used to 
fix the scale of the coulombic contributions. For the remaining short-distance 
perturbative corrections the BLM prescription |]26| was used to justify the choice 



/i = e~^^/^^Mb. Given that, the quoted uncertainty in the final as determination 
looks rather unrealistic. 

In order to make a more reliable analysis one needs to know the size of the 
higher-order perturbative corrections. Fortunately, the 0{a1) contributions to 
the correlator n(g^) have been studied recently p7|- ^] . Although a complete 
analytical calculation of these corrections is still not available, the present infor- 
mation is good enough to perform an accurate analysis of the moments A4n- 

In this paper, we present a detailed study of the relevant moments, using all 
the information on n(g^) that we are aware of. From the present experimen- 
tal data we determine the numerical values of the bottom quark mass and the 
strong coupling. Moreover, we perform a thorough analysis of the associated 
uncertainties. 

The resulting values of as{Mz) and M^ are found to be less precise than 
what was claimed in ref. ||2^. Nevertheless, they still constitute rather accurate 
determinations. The value of the strong coupling constant turns out to be in 
excellent agreement with the more precise measurements obtained from the r 
hadronic width P^^5| or from Z -^ hadrons data [^. Previous claims |^, 



that low-energy determinations of a^ result in lower as{Mz) values than higher- 
energy ones are thus unfounded. On the other side, our analysis of the T system 
provides the most precise determination of the bottom quark mass today. 

The known perturbative contributions to the moments are given in Section 2 
and the Coulomb resummation is performed in Section 3. The non-perturbative 
corrections are discussed in Section 4. Section 5 contains the phenomenological 
parameterization, extracted from the present data. The numerical analysis is 
presented in Sections 6 and 7, which use the pole mass and modified m^inmal 
subtraction schemes, respectively. A short conclusion is finally given in Section 8. 



2 Perturbation theory 



In perturbation theory the vacuum polarization n(s) can be expanded in powers 
of the strong couphng constant, 



w\s) = n(°)(s) + an«(s) + a'n(2)(s) + ... 



(2.1^ 



with a = as/Tf- Analogously, expansions for RP^iy) and A^^* can be written. 

For the first two terms, analytic expressions are available [^,|3^. Here, we 
only give R^'^\v) and R^^\v). The corresponding formulae for 11'^'') (s) and n^^''(s) 
can for example be found in refs. |HO|,ET|. 



2 ^ ' 



(2.2) 



i?(i) 



2(l + t;2)(3-t;2) 



4Li2(p)+2Li2(-p)+ln(p)(ln(l+p) + 21n(l-p) 



4t;(3 - v'^) ( ln(l + p) + 2 ln(l - p) 



- - (1 - t;)(33 - 39^; - 17^^ + 7v^) ln(p) + ^v{h- ?,v'^) 



(2.3) 



where p = (1 — f)/(l + v) and Li2(2;) is the dilogarithmic function [Q. The 
expression for R^^^ implicitly includes a factor Cp = 4/3. 



Using the integral representation (|1.4|) for A^n, one finds the following expres- 
sions for the moments: 



A4W 



with 



A 



4 
3 



3(n + l)B(5/2,n), 

B(5/2, n) I {2n + 3) A + ^ A+i - 4n + 12 

6 2 4 6 

^ri ^ (n + 1) ~ (ra + 2) ~ (72 + 3) 



1 _ 3 _ 3 

~ 2n ~ (n + 1) ~ 2(n + 2) 



(2.4) 



(2.5) 



1 



n+l 



(n + 2)B(l/2,A;) 3 



A;B(l/2,n) 



A; (2A;-1) 



(2.6) 



and B(a;,?/) being Euler's Beta function. The first order moments AA)^' are in 
agreement with the result found by Generalis [^1. 



The second-order vacuum polarization n(^)(s) is still not fully known analyt- 
ically. However, the method of Fade approximants has been recently exploited 
to calculate 11*^^^ numerically, using available results at high energies (s -^ — C)o), 

(2) 

analytical results for the first seven moments A^) for i = 1, . . . , 7 and the known 
threshold behaviour R^'^\v) for f ^ |[27|-p9|. Following the lines of ref. [^, it 



is convenient to split II^^^ according to the colour factors, 

n(2) = clnf + CACFTi%\ + CFTniiif'^ + CFTnf, (2.7) 

and to treat the four different contributions separately, because they exhibit 
different behaviour at threshold. Il\ and IItva contain the pure gluonic contri- 

f2) 

butions; the first term is already present in an abelian theory, whereas Htva i^ 

f2) (2) 

characteristic of the non-abelian aspects of QCD. The contributions n[ and Up 
arise from diagrams with internal light and heavy quark loops respectively.^ The 
spectral function R] \v) is known analytically and Rp (v) receives contributions 
from a two-particle cut with threshold at 2M which is known analytically and 
a four-particle cut with threshold at AM which can be calculated numerically 
from a two-dimensional integral |^,^. These results can be used to check the 



reliability of the Fade approximation for the moments. We shall not repeat the 

t(2) 



(2) 

technicalities of the calculation of the IIV , but refer the reader to ref. 29 for 



details. 

In table |1|, we give the first twenty moments A^^^-*, -Mn^ and A^^^-*, as well 
as the four contributions to A4^^ separately. The first seven moments for Af^n 
correspond to the analytic expressions of ref. [^, whereas the moments for n > 8 



(2) 

are our results obtained from the Fade approximants. In the case oi Ai\\^ the 
values arise from a [5/4] approximant, and the moments AtTVAn; -^in ^^^ A4pn 
were calculated from [4/4] approximants because the constant contribution to 11^^^ 
in the limit f — > is unknown. To check the stability of our results and thus the 
reliability of the Fade approximation, either different Fade approximants using 
the full set of information can be calculated, e.g. [4/5] or [6/3] in the case of 

,(2) 
'■A,ny 

'■A,n 

(2) 

moment is removed as an input datum. In particular, A1^2o changes by 0.002 

(2) (2) (2) 

and Al)vA 20 ^"7 0.0003. The moments A(; „ and Al^„ can also be calculated from 



(2) 

M.\n, or Fade approximants with one order less can be constructed by removing 

(2) (2) 

one datum. For Al\„ and AfjvAn ^^e largest change is found if the seventh 



(2) (2) 11 rn I h (2) 

the available results for Ri and Rp |30|,|31[|. In the case of Af^^^, for n > 



^With respect to the bottom quark, we consider the up, down, strange and charm quarks 
to be massless. 



n 


A^(f) 


-M« 


M'^ 


MZ,n 


A<!S 


mZ 


A^(2) 


1 


2.4000 


12.1481 


11.4197 


15.9696 


-5.2627 


1.6358 


71.2368 


2 


1.0286 


7.9822 


14.3850 


14.1999 


-4.8914 


0.6010 


69.7300 


3 


0.6095 


6.0448 


15.0503 


12.1448 


-4.2652 


0.3373 


64.1861 


4 


0.4156 


4.8899 


15.0403 


10.5729 


-3.7595 


0.2238 


59.1538 


5 


0.3069 


4.1158 


14.7923 


9.3710 


-3.3623 


0.1627 


54.9238 


6 


0.2387 


3.5585 


14.4586 


8.4283 


-3.0453 


0.1252 


51.3801 


7 


0.1926 


3.1370 


14.1001 


7.6699 


-2.7871 


0.1003 


48.3813 


8 


0.1596 


2.8067 


13.7427 


7.0464 


-2.5727 


0.0827 


45.8120 


9 


0.1351 


2.5406 


13.3977 


6.5244 


-2.3917 


0.0698 


43.5843 


10 


0.1163 


2.3214 


13.0696 


6.0804 


-2.2368 


0.0599 


41.6318 


11 


0.1015 


2.1377 


12.7599 


5.6980 


-2.1025 


0.0521 


39.9043 


12 


0.0896 


1.9815 


12.4685 


5.3648 


-1.9850 


0.0459 


38.3628 


13 


0.0799 


1.8468 


12.1945 


5.0718 


-1.8811 


0.0409 


36.9772 


14 


0.0718 


1.7296 


11.9369 


4.8118 


-1.7886 


0.0367 


35.7233 


15 


0.0650 


1.6267 


11.6945 


4.5796 


-1.7057 


0.0331 


34.5821 


16 


0.0592 


1.5354 


11.4661 


4.3706 


-1.6309 


0.0301 


33.5379 


17 


0.0542 


1.4541 


11.2505 


4.1816 


-1.5629 


0.0276 


32.5780 


18 


0.0499 


1.3810 


11.0468 


4.0097 


-1.5010 


0.0254 


31.6919 


19 


0.0461 


1.3150 


10.8539 


3.8527 


-1.4443 


0.0234 


30.8707 


20 


0.0428 


1.2552 


10.6709 


3.7086 


-1.3921 


0.0217 


30.1070 



Table 1: One-, two- and three- loop perturbative contributions for the first twenty 
moments Ain- 



the moments we are interested in, the contribution of the four-particle cut can 
be neglected, and it is sufficient to consider the analytically available expressions 

(2) 

for the two-particle cut. The agreement with the twentieth moments M.12Q and 
M.p20 ^s calculated from the Pade approximants is better than 10~^ in both 
cases. Thus for all moments under consideration the uncertainty is below 0.02%, 
being completely negligible for our application. We conclude that the method of 
Pade approximation works sufficiently well to predict the moments up to at least 
n = 20. 

Generally, the moments A^jf-* depend on the renormalization scheme and scale 



for the strong coupling constant. The values presented in table correspond to 
a{M) in the MS scheme and a renormalization scale //„ = M. If the scale fia is 
varied, the moments change according to 

-A^i'H/^a) = -Mi^'(M)+/5iA<«ln^. (2.8) 

The first coefficients of the /3-function in our conventions are given in appendix A. 

Besides varying the scale /Iq at which the coupling constant is evaluated, in 

our numerical analysis we shall also use a different definition of the quark mass. 

Apart from the pole mass, the sum rules will also be analyzed in terms of a 



running MS mass m(yUm), evaluated at a scale fim- Of course, physical quantities 
should remain unchanged. Therefore, the variation in our results originating from 
changes of scheme and scale in the coupling and quark mass will give an estimate 
of the uncertainty due to higher orders in perturbation theory. It should already 
be remarked that we have deliberately chosen different scales /i^ and fim in the 
coupling and mass respectively, in order to be able to vary them independently. 
From the definition ( |1.3[ ), it is easy to calculate the relations between the 
moments defined in terms of the pole mass and those expressed in terms of a 



running MS mass m{^m)- 

ATi'^ = A^« + 2nr«A<(f), (2.9) 

:M;f = ^(f) + 2nr«A<«+n(2ri^) + (2n-l)r«^)ML°), (2.10) 
where r^^ and r^-* appear in the relation between pole and running MS mass 

m{iXm) = M [l + a{ij,a) r^r^\nm) + a{fXaf r'^\iJa, fJ-m) + ■ ■ ■ \ ■ (2.11) 

Explicit expressions for r^^ and r^^ are also given in appendix A. 

As can be seen from table |T], for large n the higher-order corrections grow with 
respect to the leading order. At n = 8 the first order correction is roughly 120% 
of the leading term whereas the second order contribution is 140%. At ra = 20 
the contributions of first and second order are 200% and 340% respectively. This 
behaviour of the perturbation series for large moments is well known [p!7| ,p!8|,p5 



and originates from the fact that the relevant parameter in the Coulomb system 
is ag/v which leads to a dg^Jn dependence of the moments. Thus for higher n 
the perturbative corrections become increasingly more important and have to be 
summed up explicitly in order for the theoretical expressions to make sense. This 
Coulomb resummation will be discussed in the next section. 



If, on the other hand, the MS mass is used, it is not clear how a Coulomb 
resummation could be performed, because now the velocity v depends on the 



renormalization scale. However, the radiative corrections in the MS scheme are 
somewhat smaller than if a pole mass is used. At ^a = f^m = ""^ and n = 8 the 
first and second order corrections are —27% and 3% whereas for n = 20 they 
are —170% and 115% respectively. This suggests to try to find a scale /!„ for 
which the perturbative corrections stay within a reasonable range. If we require 
the second-order correction not to exceed 50%, the scale fim should lie within 
2.7 GeV ^ Hm ^ 3.7 GeV. For the numerical analysis it is therefore possible to 



also exploit the sum rules in the MS scheme if /i^ ~ 3.2 ± 0.5 GeV is chosen. 



3 Coulomb resummation 



Let us first state our Ansatz for the Coloumb resummed spectral function R{v) 
and then discuss the different components: 



R{v) 



1 - ACpa + 16C|a2) { i?W + Re + R^^^a + R^'^^a^ } , (3.1^ 



with 



Re -- 


9 
2 


Xy 

(1 - e-^^/^) 


— V 


1 




i?(^) = 


= i?«+4C^i?(°)- 






i?(2) = 


= R 


2)+4C^i?«- 


37r 




9 2^ 



(1) 



(3.2) 

(3.3) 
(3.4) 



Here, xy = Ti'^Cpav and ay is the effective coupling which corresponds to the 
heavy quark-antiquark potential. Expressed in terms of the MS coupling, we 
have 



ay{q^) = a{fia) f 1 + a{fia) r\^\q^ / fJ'l) + a{^^af r'~y\f / ^il) + 



(3.5) 



Because ay is related to the static QCD potential it is independent of the renor- 
malization scale but it does depend on the three- momentum transfer between 
the heavy quark and antiquark. Explicit expressions for Vy and Vy are given in 
appendix B. 

The term Rq corresponds to the resummed spectral function resulting from 
the imaginary part of the Green function for the QCD Coloumb potential. It 



8 



resums the leading (a/v)"' and some of the sub-leading corrections [Q. The cor- 
responding terms have to be subtracted from R^^^ and R^'^\ Although the QCD 
corrections to R are only known to order a^, we have included the recently calcu- 
lated O(a^) contribution for ay |^,0, in order to investigate the dependence of 
our results on higher-order corrections. This will be discussed further in section 6. 
In addition, we have factored out the correction to the vector current which orig- 
inates from transversal, hard gluons. To the known correction "—4(7^0" we have 
added a term IGCpa"^ in order not to generate additional corrections of order a^ 
proportional to R^'^\ We shall comment further on this point below. After per- 
forming the Coulomb resummation, the large-moment behaviour of the remaining 
terms is much weaker. Let us discuss the different contributions in more detail. 
It should be clear from eq. ( |1.4| ) that the large-n behaviour can always be 
inferred from the small-f behaviour of R{v). Expanding eq. ( |2.4| ), we obtain 

90F [ 7 145 1645 



All°) 



4^3/2 



2% 27^2 



)10^3 



+ 



n'* 



(3.6) 



The small-f expansion of R^^^ is given by: 
i?(i) = Stt^ - 24f + 2tt^v^ + ( 16 ln(8t;2^ 



(iQ ln(8t;2) - —)v^ - ttV + 0{v^) . (3.7) 



The first two terms are canceled by the additional contributions to R^^\ There- 
fore, although Ai^^ /A4'^^ increases as y/n, now 



-M« 



in 



3/2 



M 



(0) 



9v^ 



16 
3n 



1 / ^\ 11 



TT' 



3/2 ^ ^ 

+ 0\ 



n 



3/2 



n^ 



(3.8) 



For the moments n = 8, . . . , 20 the second term is of the same size as the first. 
Thus for the case of interest the large-n expansion is very badly behaved. In 
fact, below n ^ 100 the ratio tVI^^VTW^^-* increases and only for n > 100 the 
asymptotic l/y/n decrease is approached. 

(2) (2) 

The available analytical results for R] and Rp allow to calculate the small-t; 
behaviour for these functions as well: 



R 



(2) 



37r2 r Av^ 
—— m- 



11 n 

H V ^ 

2 2 



2 r 4^2 

In 



'1 



v^ 



17 

y 



R 



(2) 



An' 



9\ /245 



v^ + 0{v^ 



v^ + Oiv^) ,{3.9) 



(3.10) 



Again, the first term in eq. ( p.9|) is canceled by the corresponding piece in the 
last term of eq. (|3.4|) if we substitute q^ = v'^s = Av'^M'^ /{I — v"^) and if the first 



9 



coefficient of the /5-function in r-y is evaluated with ni hght quark flavours. On 

(2) 

the other hand, Rp vanishes at threshold and hence has no contribution which 
should be resummed in the Coulomb term. This indicates that consistently the 
coupling constant in Re should be evaluated in an effective theory with only rii 
active flavours. To facilitate the numerical analysis, we then prefer to rewrite the 
full expression for A^„ in terms of the coupling a deflned in the n^-flavour theory. 
From the matching relations for a |^5|-^, it follows that this just amounts to 
using the corresponding Pi with rii flavours in eq. (|2.8| ). 

Analogously, all terms of 0{l/v), (9(lnt>^) and C(l) for R)^ and Rj^a ^^^ 
canceled in eq. ( |3.4| ), such that w"^' vanishes in the limit f — > 0. Nevertheless, 
for these two functions the contributions of 0{v) which determine the constant 
terms in A^„^/A^.[f^ and A^„Xryi/-^n^ ^^^ not known analytically. Precisely 
those terms correspond to the current correction from transversal, hard gluons. 
In order to obtain information on the large-n behaviour of the second-order mo- 
ments, we can assume an expansion analogous to eq. (|3.8|), however including 
a constant term and fltting this Ansatz to the moments as calculated from the 
Pade approximation. We then flnd 

M^) 1174.3 Inn 534.6 

— ^ ^ 161.1 ^ + 819.9 + + . . . . 3.11 



Mn' V^ n n 

The flt has been determined using moments with n = 20, . . . , 50, but the coeffi- 
cients are rather stable if the number of flt points is changed. Although the error 
on the coefficients probably is substantial, it nevertheless shows that again here 
the large-n expansion converges slowly and for the range of n in which we are 
interested, higher-order terms have to be included. 

From the constant term in eq. ( p.ll| ), we can in principle infer the short 
distance correction resulting from transversal gluons. However, in the region of 
interest, namely for n = 8, . . . , 20, the ratio Ai^^y-^n^ ~ 40. This contribution 
should be added to the l6Cp already factorized in eq. (^]l]), therefore further 
increasing this positive correction. In the work by Voloshin |2^, the BLM scale 
setting prescription p6[ was applied to absorb the O(a^) correction in the term 
—ACp a{fia) by changing the renormalization scale /i^, and it was found that this 
should be accomplished with the choice /!„ ~ 0.63M(,. Because the flrst and 
second order terms appear with different signs, from the explicit calculation we 
now see that on the contrary the scale fia should be greater than M^. Because we 
keep the 0(0^) correction explicitly, there is no need to evaluate this scale here. 

10 



4 Gluon condensate 



Analytical results for the gluon condensate contribution to the massive vector 
correlator are available at the next-to-leading order [§,|4^. Adopting the notation 
of ref . ^^ , the corresponding moments are given by 



M 



Sn' 



n,G^ 



{aGG) y 



M4 



l+aK 



with 



,y 



-- (n + l)(n + 3)5(1/2, n + 3), 



(4.1) 



(4.2) 



and the coefficients h^ together with numerical values for the ajf up to n = 21 
are shown in table ^ The coefficients 6^ depend on the renormalization scheme 
for the mass. If the MS scheme is used the h^ change according to 



K 



hi + (2n + 4) r. 



(1) 

m 



(4.3) 



n 


1 


2 


3 


4 


5 


6 


7 


< 


-0.3048 
10.4768 


-0.5079 
11.7202 


-0.7388 
12.8494 


-0.9946 
13.8928 


-1.2730 
14.8685 


-1.5726 
15.7888 


-1.8918 
16.6625 


n 


8 


9 


10 


11 


12 


13 


14 


< 
hi 


-2.2296 
17.4964 


-2.5851 
18.2956 


-2.9574 
19.0643 


-3.3457 
19.8058 


-3.7495 
20.5229 


-4.1682 
21.2179 


-4.6012 
21.8928 


n 


15 


16 


17 


18 


19 


20 


21 


bl 


-5.0482 
22.5492 


-5.5087 
23.1887 


-5.9823 
23.8124 


-6.4686 
24.4216 


-6.9674 
25.0173 


-7.4784 
25.6002 


-8.0012 
26.1712 



Table 2: First- and second-order coefficients for the gluon condensate contribu- 
tion to the moments Ain- 

From eq. ( |4.2| ) it is clear that the relative growth of A4^Q2/Ain^ is propor- 
tional to n^. Therefore, the non-perturbative contribution grows much faster 
than the perturbative moments. In addition, as can be seen from table ^ in the 
pole-mass scheme at yUa = M the next-to-leading order correction is of the same 
size or larger as the leading term. Because the perturbative expansion for the 
gluon condensate cannot be trusted, we shall restrict our analysis to a range of 
n where its contribution to the moments is small and can be neglected. Using 
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{aGG) ^ 0.021 GeV^ HIH, we find that for n < 20 tlie contribution from the 
gluon condensate to the hb moments is below 3%. Thus, we shall restrict our 
phenomenological analysis to this range. 



In the MS scheme the situation concerning the perturbative expansion is 
somewhat better. If we take /im ~ 3.2 GeV, as was discussed at the end of 
section 2, for n < 20 the next-to-leading order contribution stays below 70% of the 
leading order. This demonstrates that a determination of the gluon condensate 
from charmonium should be performed in the MS scheme. Nevertheless, for this 
work also in the MS scheme we shall keep the restriction to n < 20. A lower 
limit on the number of moments will be discussed in the next section. 



5 Phenomenological parameterization 

In the preceding sections, theoretical predictions for the spectral function R[s) 
and the related moments M.n have been calculated without further specifying the 
actual quark content. For the phenomenological parameterization of the spectral 
function we shall now restrict our discussion to the hb system. 

In the narrow-width approximation the contribution to Rb{s) from a T{kS) 
resonance is given by 

RbMs) = ^r{T{kS)^e+e-)MksS{s-M^,s), (5-1) 

where a denotes the running QED coupling evaluated at a scale around the 



resonance mass. Because in the Review of Particle Properties |^0[ the electronic 
widths have been calculated with a^ = 1.07 a^ where a = 1/137.04 is the fine 
structure constant, we shall use this value accordingly. In the case at hand the 
narrow-width approximation is extremely good because the full widths of the 
first three T resonances are roughly a factor 10"^ smaller than the corresponding 
masses and the higher-resonance contributions to the moments are suppressed. 

Experimentally, the first six resonances have been observed. The measured 
masses and electronic widths are collected in table 0. For our numerical analysis 
the errors on the masses can be safely neglected and have thus not been listed. 
Inserting eq. ( p.l[ ) in the definition of the moments A4n, eq. ( |1.3| ), we obtain 

M.-,M,r{^t^.Ps!^}. ,5.) 

So 
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k 

Mks [GeV] 
Tks [keV] 



12 3 4 5 6 

9.460 10.023 10.355 10.580 10.865 11.019 

1.31 ±0.04 0.52 ±0.03 0.48 ± 0.08 0.25 ± 0.03 0.31 ± 0.07 0.13 ±0.03 



Table 3: Masses and electronic widths of the first six T{kS) resonances. 



The numerical weight of the heavier resonances in ( |5.2| ) decreases strongly for 
increasing values of n. The contribution of the T(5S') [T(65')] state is 9.5% [4%] 
at n = 0; 1% [0.3%] at n = 10; and a tiny 0.08% [0.02%] at n = 20. Therefore, 
taking n ^ 10, the uncertainties associated with the contributions of higher-mass 
states are very small. 

The second term in eq. (|5.2| ) accounts for the contributions to Rh above the 
sixth resonance and is approximated by the perturbative continuum. Generally, 
the continuum threshold ^/sq should lie around the mass of the next resonance, 
which has been estimated in potential models |^. For our analysis we shall use 



y/so = 11.2 ± 0.2 GeV. The lower value for sq includes the mass of the sixth 
resonance and should be a conservative estimate. There is still a contribution 
missing which stems from open B production above the BB threshold and below 
Sq. From the experimental data [^ we infer that its infiuence is small and has 



been included in the variation of sq. 

6 Numerical analysis in the pole-mass scheme 

Quark-hadron duality entails the equality of the theoretical moments Ai*^ pre- 
sented in sections 2 to 4 and the phenomenological moments Ai^^ discussed in 
the previous section. The moments corresponding to the Coulomb term Rq of 
eq. ( |3.2| ) have been calculated from eq. ( |1.4| ) by numerical integration. To sup- 



press higher resonances as well as power corrections, following ref. [^, we have 
restricted n to the range n = 8, . . . , 20. Solving the moment sum rules for M^, we 
can fit Mb to a constant by varying Mb and as{Mb). The fit has been performed 
using the program Minuit ||5^ . For the central set of parameter values our result 
is 

Mb = 4.604 ±0.009 GeV, (6.1) 

a,(Mfe) = 0.2197 ±0.0097. (6.2) 
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The error in these results just corresponds to the statistical error of the fit. In 
the fit we have included every second moment to have less statistical dependence, 
but the results change very little if all moments with n = 8, . . . , 20 or only every 
fourth moment is used. In figure ^ the resulting values for Mt are displayed as a 
function of n. This illustrates that a constant Mf, in the range 8 < ra < 20 really 
produces an excellent fit. 




Figure 1: The b quark pole mass as a function of n. 



In the remaining part of this section let us present a detailed discussion of the 
errors resulting from the various input quantities. A compilation of all different 
contributions to the errors on Mf, and as{Mf,) is summarized in table ^ The 
dominant theoretical uncertainty is due to the unknown higher-order perturbative 
corrections. We have estimated this uncertainty in three different ways. As has 
been already remarked in section 3 the relation between the effective coupling 
ay in the Coulomb potential and a*^'^ is known to 0{a^) ||4^,Q. We can thus 
include this correction in Re to see what the influence on our results is. Although 
this is not consistent because the corresponding correction R^^^ is not available, 
it nevertheless can be taken as an error estimate of higher-order corrections. A 
second possibility is a variation of the scale at which a^ is evaluated. For the 
result in table ^ we have chosen the range Mh/2 < fia ^ 2Mh. As a final test on 
the importance of higher-order corrections, we can remove the 0{a'^) term i?^^^ 
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AMb [MeV 


Aas{M,) [10-3] 


statistical 


±9.1 


±9.7 


0{a^) Coulomb 


±7.3 


±22.2 


0{a') 


±0.9 


±8.2 


scale fia 


i2.1 
=^1.3 


_|_28.7 
=■=19.6 


continuum 


i2.6 
=■=2.0 


1 3.7 
="=2.7 


{aGG) 


±5.3 


±3.6 


Te+e- 


±3.1 


±6.7 


total 


±13.5 


±26.9 



Table 4: Separate contributions to the errors of Mf, and as{Mh 



completely. From table ^ we observe that including the 0{a^) correction to ay 
has a much bigger influence than removing w'^^ . This is not unexpected because 
the Coulomb piece sums up the dominant contributions in the large-n limit. The 
uncertainty of the scale dependence is of the same order as the sum of the other 
two contributions. For our estimate of the uncertainty resulting from higher- 
order corrections we can now either take the scale dependence or combine the 
other two contributions. Adding all three would double count the error, because 
the uncertainty in an asymptotic series, such as the perturbative expansion, is 
bounded by the size of the last known term. For our flnal results, we have chosen 
to include the errors of varying the Coulomb and the O(a^) terms. 

The error from the continuum contribution has been estimated by varying Sq 
in the range ^/sq = 11.2 ±0.2 GeV. The entry for the gluon condensate in table ^ 
results from removing the gluon condensate completely and for the uncertainty 
from the electronic widths we have varied all widths within the errors given in 
table 1^. With respect to the uncertainty resulting from higher orders all these 
errors are small. Adding all errors in quadrature, we arrive at our flnal result in 
the pole-mass scheme: 



Mb = 4.604 ±0.014 GeV; 
asiMb) = 0.2197 ±0.0269. 

Evolving the strong coupling constant to Mz, we flnd 



(6.3) 
(6.4) 



as{Mz 



0.1184 ± 



0.0070 
0.0080 



(6.5) 
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Our central result is in astonishingly good agreement to the current world av- 



erage [0, although the error turns out to be larger. Further comments on our 
results also with respect to the paper by Voloshin [^] have been relegated to the 
conclusions. 



7 Numerical analysis in the MS scheme 

Besides analyzing the moment sum rules exploiting the pole mass M^, in addition 



we have investigated the same sum rules in the MS scheme. In contrast to the 



pole mass, the quark mass in the MS scheme depends on the renormalization 
scale /im- As has been remarked in section 2, to restrict the O{o?') corrections to a 
reasonable size, Hm should lie in the range /Xm = 3.2±0.5 GeV. We have refrained 
from performing a resummation of the large radiative corrections because now the 
velocity v depends on the renormalization scheme and it is not straightforwardly 
possible to proceed in analogy to the Coulomb resummation for the pole mass. 

The fitting procedure was performed along the same lines as for the pole-mass 
case. For the central values of our input parameters, we obtain 



nxbinik) = 4.133 ±0.002 GeV, 
asirrib) = 0.2325 ± 0.0044 , 



(7.1) 
(7.2) 



where again the errors are purely statistical. Since it is more standard to eval- 
uate the running 6-quark mass at rrif,, we have evolved our immediate result 
mh{3.2 GeV) to this scale with the help of the renormalization group equation. 





Arrib [MeV 


Aa,(mb) [10-3] 


statistical 


±2 


±4.4 


scale Hm 


_i_33 
=■=36 


_i_36.2 
=•= 5.8 


scale fia 


±^? 


1 25.1 
=■=28.0 


continuum 


±1 


1 3.5 
=■=2.6 


(aGG) 


±2 


±2.3 


Te+e- 


±3 


±6.0 


total 


_|_59 


_|_44.9 
=■=29.7 



Table 5: Separate contributions to the errors of mb{mb) and asirrib) 
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The separate contributions to the theoretical error have been obtained by 
performing the same variations as for the pole-mass scheme and have been listed 
in table ^. The uncertainty from higher-order corrections is now due to the 
variation of the scales /im = 3.2 ± 0.5 GeV and 2.6 GeV < fia ^ 2mb. The scale 
fia should not be taken lower than roughly 2.6 GeV because otherwise the O(a^) 

(2) 

correction Ai^ becomes unacceptably large. Since in addition to j2a for the MS 
scheme we can also vary /x^, the resulting uncertainty, especially for m^, is larger 
than for the pole mass. Adding all errors in quadrature, we arrive at our final 
result in the MS scheme: 

rribimb) = 4.13 ± 0.06 GeV, (7.3) 

asimb) = 0.2325 ±°:°^f7 . (7.4) 

Evolving the strong coupling constant to Mz, we find 

asiMz) = 0.1196 ±Ho80- (7-5) 

It is gratifying to observe that the resulting values for as{Mz) from the pole- 
mass and MS schemes turn out to be in very good agreement. This is a further 
indication that the uncertainty from unknown higher-order corrections is under 
control. In addition, our results mb{mb) and Mb for the 6-quark mass satisfy the 



relation ([A.4|) between the pole and MS mass within the errors. This should be 



expected because the relation ( |A.4D has been used to rewrite the moment sum 



rules in terms of the MS mass. Nevertheless, it again shows that variations due 
to higher orders are accounted for by our error estimates. 

8 Conclusions 

Before we enter a discussion of our findings, let us again summarize the central 
results. For the bottom quark mass in the pole-mass as well as MS scheme, we 
obtain 

Mb = 4.60 ±0.02 GeV, (8.1) 

mbimb) = 4.13 ±0.06 GeV, (8.2) 

respectively. Combining both determinations of the strong coupling constant ««, 
we find 

as{Mz) = 0.119 ±0.008. (8.3) 
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We have not averaged the errors of the two determinations because they are not 
independent. 

The bottom quark mass values obtained by us are in good agreement to 



previous determinations from QCD sum rules p2|, |23|, |5^-p6| and a very recent 
calculation from lattice QCD ||5^. Owing to the big sensitivity of the moment 
sum rules for the T system to the quark mass, and the good control over higher- 
order Us corrections, our result is more precise. 

Nevertheless, the pole quark-mass value obtained by us is in disagreement to 
the result found by Voloshin p4|. In our opinion the discrepancy is due to the 



importance of higher 0{l/n) corrections, which in ref. |2^ were either neglected 



or numerically fitted from the sum rules. In |24] it was assumed that the leading 



order correction goes like 1/n. However, from eqs. ( |3.8| ) and ( |3.11 ), it is clear 



that they rather behave like Xj \fn. Besides, we have also demonstrated that for 
the region of n used in the analysis, the large-n expansion is not justified. In 
addition, the second-order a^ correction was only partially and partly incorrectly 
taken into account. Therefore, the scale dependence of a^ was not under control. 
Let us shortly comment on the renormalon ambiguity of the pole mass. Dur- 
ing the last years, it has been realized that beyond perturbation theory the 
pole masses for the charm and bottom quarks suffer from unknown renormalon 
ambiguities, leading to additional theoretical uncertainties in their determina- 
tion |S^-^. On general grounds this uncertainty has been estimated to be of 
(9(100 Me V). Throughout our analysis, the pole mass has been defined as the 
pole of the perturbatively renormalized quark propagator. Our determination 
(|8.1|) might therefore be subject to additional uncertainties which go beyond per- 
turbation theory but which we cannot assess in a precise way. 



Within our errors the result obtained for ^^(M^), eq. (|8.3| ), is compatible 



with the result by Voloshin P3 , though, given the shortcomings of this analysis 



discussed above, his errors appear to be largely underestimated. Thus, previous 
claims of a low value of agiMz) from low-energy determinations which could hint 
to new physics P3,P^ are unfounded. On the other hand, our central value for 



agiMz) is surprisingly close to the current world average as{Mz) = 0.118 ±0.003 
[^, although the error is certainly larger. 

The dominant uncertainty for the determination of the fe-quark mass and a^ 
from the T system was found to originate from the dependence on the renor- 
malization scale, or, equivalently, the size of the as yet unknown higher-order 
corrections. Improving the error on the a^ determination will thus only be possi- 
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ble if the full 0{al) correction to the moments Ain is known and if it turns out 
to be reasonably small. We hope to return to this question in the future. 
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Appendices 



A Renormalization group functions 

For the definition of tlie renormalization group functions we follow the notation 
of Pascual and Tarrach ||6l| , except that we define the /3-function such that (3i is 
positive. The expansions of /9(a) and 7(a) take the form: 

/3{a) = - Pia-P^a^-Psa^-- ■ ■ , and 7(a) = 710+720^+730^ + ..., (A.l) 



with 

A - 
and 



1 r 



HCa-^Tti 



7i 



C, 



132 = ^[nCl-lQCATuf-QCFTuf], (A.2) 



72 



48 



97Ca + 9Cf - 20Tr2 



The relation between pole and running MS mass is given by 



m{^r 



M 



1 + a{^a) r^m (/^m) + 0■{^^af T^m (/^a, ^^m) + 



where 



.(1) 



.(2) 



^m?o - 7i In 



,(2) 
m,0 



72 + (71 - Pi) ^m?0 



(A.3) 



(A.4) 



(A.5) 



In^ + — (71 - /5i) In —— 



7i + iJiin — 

fJ-a 



.(1) 



(A.6) 



The coefficients of the logarithms can be calculated from the renormalization 
group and the constant coefficients rl^Q and rl^Q are found to be p9|,|6^,|63[| 



,(1) 

m,0 



„(2) 
' m,0 



-c. 



(A.7) 



^Ki^ - f ^^^) - i^(^) + ^^(2) '-') + ^-^^KS + \^^'^) 



+ CaCA -^ + ^C(2) + ^C(3) - ^C(2)ln2) +C^t(^ - ^C(2) ) . (A.^ 
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B The effective coupling ay 



In terms of the MS coupling the effective couphng ay is given by 

ay(g2) = a{fia) [l + a{fia) r^v i^^ / f^l) + (^if^a)^ rv\^^ / f^l) + ■ ■ ■ 
where 



V 



,(2) 
V 



(1) A, r 

^ fa 



^(2) _ 



— + Pl'^V.O 



/3? --' 



1 ? , Pi 1 2 9 

1^ — + ^1^ :;¥ 



/^a 



/^a 



(B.l) 

(B.2) 

(B.3) 



Like in eq. (|A.4| ) the coefficients of the logarithms are determined by the renor- 



malization group and the constant coefficients ryj, and ryl are found to be 



[ggigg 



„(i) 

Vfl 



.(2) 
V,0 



C /I I Til . 

36 9 

1 [^2/4343 
16r^ll62" 



(B.4) 



r 2 ^\ 22 \ ^1798 56 A 



- CpTni 



55 



16C(3) 



T^nf 



400 



3 ''V ' 8f J 

Here, n^ = nj — 1 is the number of light quark flavours. 



(B.5) 
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